Dynamical properties of hard-core anyons in one-dimensional optical lattices 
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We investigate the dynamical properties of anyons confined in one-dimensional optical lattice 
combined with a weak harmonic trap using the exact numerical method based on a generalized 
Jordan- Wigner transformation. The density profiles, momentum distribution, occupation distribu- 
tion and occupations of the lowest natural orbital are obtained for different statistical parameters. 
The density profiles of anyons display the same behaviors irrespective of statistical parameter in 
the full evolving period. While the behaviors dependent on statistical property are shown in the 
momentum distributions and occupations of natural orbitals. 
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I. INTRODUCTION 

Generally according to the quantum statistical prop- 
erty particle is classified as boson and fermion. The 
wavefunction of identical bosons is symmetric under ex- 
change while that of identical fermions is antisymmetric 
under exchange. As a natural generalization, physicist 
proposed that there exists anyon interpolating between 
Bosc and Fermi statistics, which satisfies fractional statis- 
tics [l|. It has become an important concept in con- 
densed matter physics [HH and has ever been used for 
successfully explaining the fractional quantum Hall ef- 
fect (FQHE) [5(. Fractional statistics also play impor- 
tant roles in the theory for non-fermi liquid, Chern-Simon 
theory, and other aspects @, 0] ■ Another potential ap- 
plication of anyons is to realize the topological quantum 
computation with non-abelian anyons [H, Q . Besides the 
traditional two dimensional electron system, cold atoms 
in low dimension is also a popular platform to realize the 
fractional statistics. It has been suggested that anyons 
can be created, detected and manipulated in rotating 
Bose-Einstein condensates (BECs) and cold atoms in op- 
tical lattices [IM1- Particularly, researches on anyons 
in cold atoms are not restricted in two dimensional sys- 
tem but the proposal to realize anyons in one dimensional 
(ID) optical lattice is also put forward recently [l3j |. 

Since BECs are realized experimentally, cold atoms has 
been making great progress both experimentally and the- 
oretically for its "purity" and high controllability com- 
paring with the traditional condensed matter system. 
For the dimensionality effect and extremely profound 
strong correlation in low dimensional system, the quan- 
tum gas in low dimension has also attract much at- 
tention [l4j |. With anisotropic magnetic trap or two- 
dimensional optical lattice potentials, the particle mo- 
tion is tightly confined in two directions to zero point 
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oscillations [TH-ITtJ and the strongly correlated Tonks- 
Girardeau (TG) regime can be achieved [HI, fl7j | ■ By 
crossing the confinement-induced resonance (CIR) from 
the TG gas, the super TG (sTG) gas is accessible, which 
show stronger correlation than TG gas [l8j|. The in- 
teraction between atoms can be tuned in the full in- 
teracting regime with Feshbach resonance technique and 
confinement-induced resonance by tuning magnetic field. 
It is the excellent tunability of cold atom that makes it to 
be an excellent candidate to realize fractional statistics. 

In factjbefore the experimental scheme of Keilmann 
et. al. [13|, ID anyon gas has been investigated the- 
oretically in various ID systems ITT? 22 including the 
Bose quantum gas. Kundu proved that a ID Bose gas 
interacting through ^-function potential combined with 
double (5-function potential and derivative <5-function po- 
tential is equivalent to the anyon gas interacting via 5- 
function potential [2l|. This stimulated many research 
interests on (5-anyon gas [22- 32] . It turns out that the 
ground state density distribution of 5-anyon gas displays 
similar behavior as that of Bose gas with the increasing 
interaction. In the strong interacting regime the density 
distribution shows the same behavior as the free fermion 
[33M37| . which is irrespective to the statistical parame- 
ter. The special property resulted from the fractional 
statistics exhibits in the reduced one body density ma- 
trices and the momentum distributions [38l . |39| . The 
momentum distribution of anyon differs from fermion's 
oscillations and boson's single peak structure, which are 
symmetric about the zero momentum. The momentum 
distribution of anyons is asymmetric when the statisti- 
cal parameters deviate from the Bose and Fermi limit 
[27H3l| . This special behavior originates from that the 
reduced one body density matrix is a complex Hermitian 
one rather than a real one for Boson and Fermion. 

While most studies have focused on the static proper- 
ties of the ID anyonic gas, its dynamics remains to be 
investigated. The present paper shall study the dynam- 
ics of anyons confined in optical lattice with weak har- 
monic trap in the hard core limit. In Ref. [28j . we have 
extended the exact numerical method originally used to 
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treat hard-core bosons by Rigol et al. [40( to deal with 
the ground state of hard-core anyons (HCA) in optical 
lattice. Here, we further extend this method to inves- 
tigate the dynamics of HCA. By evaluating the exact 
time-dependent one-particle Green's function, we obtain 
the reduced one body density matrix (ROBDM) and thus 
the density profiles and the momentum distribution at 
arbitrary time. The dynamical properties induced by 
anyonic statistics shall be displayed in the momentum 
distribution. 

The paper is organized as follows. In Sec. II, we give 
a brief review of ID anyonic model and introduce the 
numerical method. In Sec. Ill, we present the density 
profiles, momentum distributions, occupation distribu- 
tion and the occupations of the lowest natural orbital for 
different statistics parameters. A brief summary is given 
in Sec. IV. 



II. FORMULATION OF THE MODEL AND 
METHOD 

We consider N hard core anyons confined in an optical 
lattice of L sites with a weak harmonic trap and the 
second quantized Hamiltonian can be formulated as 
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Here the harmonic potential Vi = V (l- (L + l)/2) 2 with 
the strength of the harmonic trap Vq. The anyonic op- 
erator a\ (a{) create (annihilate) an anyon on site I and 
satisfies the generalized commutation relations 
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where fj (fj) is creation (annihilation) operator for 
fermions. The above Hamiltonian of system with N 
anyons is transformed into a fermionic one with Np = N 
fermions 



h f = -*E (fUifi + H - c ) + E vtfUi 
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where the Fermionic operator satisfy the Fermi anti- 
commutation relation 
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The original question about anyons now can be inves- 
tigated by solving the model on the polarized fermions 
in optical lattice. We can obtain the exact many body 
wavefunction of polarized fermions with the diagonalized 
method and therefore the ground state and interesting 
physical phenomena of hard core anyons. 

The equal-time Green's function for the hard core 
anyons at time r should be expressed as 
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for j j£ I, where the sign function e(x) gives -1, or 1 
depending on whether x is negative, zero, or positive and 
X is the parameter related with fractional statistics. The 
generalized commutation relations reduce to fermionic 
commutation for x = and reduce to Bose commutation 
for x — 1; while for anyons satisfying fractional statis- 
tics x changes in between them. The hard core interac- 
tions between anyons restricts the additional condition 
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f . In the Hamiltonian t de- 



notes the hopping between the nearest neighbour sites, 
which can be tuned by changing the strength of optical 
lattice. 

In order to solve the model of hard core anyons, we ex- 
tend the numerical method to investigate the hard core 
bosons in optical lattice developed by Rigol et. al.. We 
can map the above model into the polarized fermionic 
Hamiltonian using the generalized Jordan- Wigner trans- 
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Here I^hca (t)) is the wavefunction at time r of hard 
core anyons in a system with Hamiltonian Hhca, and 
\^f{t)} is the corresponding one for the equivalent po- 
larized fermions. For the polarized fermions the time 
evolution of their initial wavefunction 
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While the matrix representation of initial wavefunction 
can be expressed as 



iff 
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Thus the fermonic time-dependent wavefunction can be 
expressed as an L x Nf matrix P(r). The matrix P (r) 
can be evaluated as 



e -irH F /hpI _ jj e -irT) 
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where U is an unitary transformation diagonalizing the 
Hamiltonian Hp, i.e., U^HpU = D with diagonal matrix 
D. After an easy evaluation the state reads 
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with 



P£ (r) = exp Hxtt) P ln (r) for I < j - 1 
Kt (t) = Pm (r) for I > j 

for n < N f and Pj£ /+1 (r) = 1 and P$ f+1 (r) = 

(I =/=j). The state \^p) has the same form with the 
replace of j by The time-dependent Green's function 
is a determinant dependent on the L x (Nf + 1) matrixes 
P' A (r) and P' s (r) 

G^(r) = <^(r)|*f (r))=det[(P'' 4 (r)) T P' s (r) 

In the present paper we will focus on the time evolution of 
the density profile and momentum distribution for hard 
core anyons with different statistical parameter \. The 
ROBDM can be evaluated by Green's function 

Pjl (r) = (a}o«) = 5 jt (1 - G iX (r)) - (1 - S^e^Gn (r) 

The diagonal part of ROBDM is the density profile and 
its Fourier transformation is defined as momentum dis- 
tribution 
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The natural orbitals (j/ 1 are defined as the eigenfunctions 
of the one-particle density matrix 



L 
3 = 1 
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and can be understood as the effective single-particle 
states with occupations A^. 
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FIG. 1: (color online) The evolving of density distribution 
for 50 hard core anyons in optical lattice of 500 sites. Vq = 
1.0 x 10 _3 t and Vb = 1.0 x lQ~ s t. 



III. DYNAMICS OF DENSITY PROFILE AND 
MOMENTUM DISTRIBUTION 

In the present paper we investigate the dynamics of 
hard core anyons in optical lattice. Initially we con- 
fine the anyons in optical lattice with a strong harmonic 
trap and then turn off the harmonic trap or reduce the 
strength of harmonic trap. The trapped anyons will 
evolve in the optical lattice. For convenience we set lat- 
tice constant a as 1, the unit of A; is 1/a and the unit of 
time is Ti/t. 

In Fig. 1,2,3 and 4 we show the evolution of 50 anyons 
in optical lattice of 500 sites. Initially the strength of har- 
monic trap Vq = 1.0 x 10~ 3 i and at r = the harmonic 
trap is released to Vb(0) = 1.0 x 10~ 8 £, i.e., the anyons 
are trapped only in optical lattice. The evolution of den- 
sity distribution is shown in Fig. 1. It turns out that 
the density distribution of anyons shall not display any 
different dynamical property from those of bosons and 
fermions. We cannot distinguish the statistical proper- 
ties by the density distribution in real space. Similar to 
hard core bosons and polarized fermions, anyons expand 
in the optical lattice and they shall gradually populate 
in the full lattice. The density of anyons at the center 
reduces faster than the density at the border because the 
anyons locating in the sites of high density posses higher 
energy. As evolving time is long enough anyons homoge- 
neously distribute in the middle regime while its density 
distribution shows peaks at the border regime. 

The corresponding momentum distribution is dis- 
played in Fig. 2. Initially bosons and fermions distribute 
symmetrically about the zero momentum as respective 
statistical property and anyons (0 < % < 1) exhibit the 
asymmetrical momentum distribution. After the har- 
monic trap is turned off for anyons of different statistical 
parameter their momentum distribution show different 
evolving properties. Fermions (x = 0.0) do not show 
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FIG. 2: (color online) The momentum distributions for 50 
hard core anyons in lattice of 500 sites, Vq = 1.0 x 10 
V = 1.0 x 10~ 8 i. (a) t = 0, (b) t = 20, (c) r = 

T = 100. 
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FIG. 3: (color online) The occupation distributions for 50 
hard core anyons in lattice of 500 sites, Vq = 1.0 x 10 _3 r. and 
Vo = 1.0 x 10~ 8 t. (a) t = 0, (b) r = 100. 



obvious change of momentum distribution (In fact the 
slight change happens according to the numerical data). 
As the statistical parameter deviates from \ = the mo- 
mentum distribution of anyon evolves from the asymmet- 
rical structure of single peak to the structure similar to 
that of Fermions but with two asymmetrical peaks. As 
statistical parameter increases these two peaks become 
more and more obvious. For the case of x = 1 (hard 
core bosons) the momentum distribution evolve to the 
structure of symmetrical double peaks. When the evolv- 
ing time is long enough the momentum distribution of 
anyons of any statistical parameter (0 < % < 1) exhibits 
identical behavior in the middle regime. 

In Fig. 3 we show the occupation distribution of nat- 
ural orbitals for the same system as above. In the full 
evolution period the occupation distribution of each or- 
bital does not change qualitatively. In the Bose limit 
anyons occupy the lower natural orbitals and the occu- 
pation distribution displays the single peak structure. In 
the Fermi limit each anyon occupies one natural orbital 
and at any time the occupation distribution seems like a 
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FIG. 4: (color online) The occupations of the lowest natural 



orbital for 50 bosons in lattice of 500 sites. Vq 
and V = 1.0 x W~ s t. 



0.8 



1.0 x W~H 



0.6 



s~0.4 



0.2 



0.0 





-1 1 1 1 1 1 1 — 


-i « 1 « 

T=0.0 






t=20.0 ■ 






- - • -x=40.0 






t=60.0 - 






t=80.0 




x=100.0 ■ 






t=120.0 

\ x=140.0" 




, . ti ; 


s 

jv£ 




• a *' ih i i i 


1 '•' v. \ 



50 



100 



150 

/ 



200 250 



300 



FIG. 5: (color online) The density distribution of 50 hard core 
anyons in optical lattice of 300 sites. Vo = 1-0 x 10 -3 i and 
Vo = 2.0 x 10~ 4 t. 



step-function. While for anyons in between these two lim- 
its, the occupations of higher natural orbitals increase as 
the statistical parameter evolves from Bose limit to Fermi 
limit. We also display the evolution of the occupation of 
the lowest natural orbital in Fig. [4] It is shown that the 
occupation is time independent in the Fermi limit and as 
deviating from the Fermi limit the occupation increases 
during the time evolution. The increase is more obvious 
for the bigger statistical parameter. When the evolving 
time is long enough the occupation of the lowest natural 
orbital shall preserve a constant. 

If the initially confined hard core anyons in optical lat- 
tice evolve in a weak harmonic trap rather than turning 
off the harmonic potential, the situation will be different. 
Anyons with different statistical parameters always ex- 
hibit the same density distributions and we also cannot 
determine the statistical properties according to evolv- 
ing property of density profiles. In Fig. [5] we display 
the evolving density distribution of 50 anyons in optical 
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FIG. 6: (color online) The momentum distribution of 50 
anyons in optical lattice of 300 sites. Vq = 1.0 x 10 — 3 i and 
V = 2.0 x l(T 4 i. (a) t = 0, (b) r = 30, (c) r = 60, (d) 
r = 80, (e) r = 100, (f) r = 120. 



lattice of 300 sites combined with a weaker harmonic po- 
tential (Vo = 2.0 x 10~ 4 <). Initially, anyons distribute 
in the middle regime of the harmonic trap and after the 
harmonic potential becomes weak anyons shall expand 
firstly. The central density decreases faster than the 
boundary density such that the density profiles behave 
as a Fermi-like distribution at r = 60. Then anyons shall 
contract because of the confinement of harmonic trap. 
With the time evolution anyons redistribute in the mid- 
dle regime of the trap. 

During the expansion, momentum distributions show 
rich dynamical structure (Fig. [6]). In the Bose limit, 
the momentum distribution firstly evolves from the orig- 
inal structure of a single peak to the structure of double 
peaks, and then back to the structure of a single peak at 
t = 60. During the later period momentum distribution 
shall behave as the single peak and double peaks alter- 
nately. In the Fermi limit the momentum distribution 
does not keep its original profile as shown in Fig. [2] and 
shall exhibit the oscillating behavior. Contrary to the 
evolution of density distribution, the momentum distri- 
bution contracts firstly and at t = 60 (Fig. [BJ:) displays 
the step- function profile in the region of — ir/4 < k < ir/4. 
Then it gradually expands to the region of higher momen- 
tum and at r = 120 (Fig. [Bf) momentum distribution 
almost recovers to the behavior at r = 0.0. For anyon 




FIG. 7: (color online) The occupation of the lowest natural 
orbital for 50 anyons in optical lattice of 300 sites. Vq = 
1.0 x 10" 3 i and V = 2.0 x 10 _4 t. 



gas in between (0 < % < 1) its asymmetric momentum 
distribution also shows the alternate behavior. For the 
case of bigger statistical parameter asymmetric double 
peaks are displayed and for the case of smaller statistical 
parameter momentum distribution exhibits the behavior 
similar to those of fermions. 

The evolution of occupation for the anyons in optical 
lattice with weak harmonic potential is similar to the case 
only confined in optical lattice. The occupation distribu- 
tion always exhibits the same behaviors qualitatively as 
those at initial time. But the occupation of the low- 
est natural evolve with time, which is displayed in Fig. 
[7] It is shown that in the Fermi limit the occupation is 
time-independent and as deviating from the Fermi limit 
it shall oscillate with the time evolution. It does not 
increase monotonously rather than fluctuates with the 
time evolution. The bigger the statistical parameter, the 
stronger the oscillation amplitude. 

In order to investigate the dynamical properties of 
anyons in Mott regime, we initially prepare a Mott state 
by superimposing optical lattice with a tight harmonic 
potential and then turn off it at r = 0.0. The dynamical 
evolutions of Mott state are displayed in Fig. [51 Ini- 
tially, the momentum distribution for anyons of any sta- 
tistical parameter exhibits the behavior similar to that 
of fermions. There arc only tiny differences in the regime 
close to zero momentum. After the harmonic trap is 
turned off, the momentum distribution of anyons in the 
Fermi limit always preserve its initial profile. Anyons de- 
viating this limit still display Fermi-like momentum dis- 
tribution in the full momentum regime, but the structure 
of double peaks appears at particular regime. In the full 
evolving period anyons in Bose limit show the sharpest 
peaks and as statistical parameter decreases (close to the 
Fermi limit) the peaks become smaller. It is same as 
before that anyons always exhibit the asymmetrical mo- 
mentum distributions except of the Bose limit and Fermi 
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FIG. 8: (color online) The momentum distribution of 50 
anyons in lattice of 300 sites. V = 1.0 x 10 -2 t and Vo = 
1.0 x 10 -8 £. (a) t = 0, (b) r = 10, (c) r = 30,(d) r = 60. 
Inset: Density distributions for the same system. 



limit. We also display the evolving density profiles in 
the inset, which is independent of the statistical parame- 
ter. At the beginning anyons are confined in the central 
regime with an anyon per site. After the harmonic trap is 
turned off, anyons expand in the lattice and peaks appear 
at the boundary regime. 



IV. CONCLUSIONS 

In summary, in the present paper we have developed 
the exact numerical method to deal with the dynamics 
of hard-core anyons confined in the optical lattice super- 
imposed with a weak harmonic potential. By evaluat- 
ing the exact time-dependent one-particle Green's func- 
tion, we obtain the density profiles, momentum distri- 
butions, occupation distribution, and the occupation of 
the lowest natural orbital in the full evolving period. It 
is shown that the evolving property of density profiles 
is independent on the statistical parameter of anyons. 



As the harmonic trap is turned off anyon shall expand 
in the optical lattice and distribute gradually in the full 
lattice. The dynamical properties of momentum distri- 
butions and occupation distributions of natural orbitals 
depend on statistical parameter of anyons. In the Bose 
limit, momentum distribution always displays the struc- 
ture of double peaks. As deviating from the Bose limit, 
the double peaks become smaller and in the Fermi limit 
the double peaks disappear. When the evolving time is 
long enough, anyons display the same momentum distri- 
bution irrespective of statistical parameter in the regime 
nearby zero momentum. If the anyon gas is initially in 
the Mott regime, during the full evolving period it shall 
be always in the Mott regime but at the particular mo- 
mentum position double peaks appear. When the har- 
monic trap is replaced with a weaker one, the density 
profiles of anyons shall exhibit the breathing behavior. 
The momentum distributions of bosons alternatively dis- 
play the structure of single peak and the structure of 
double peaks, while fermions contract and expand in the 
/c-space with the time evolution. For anyon interpolating 
between these two limit, momentum distribution also re- 
peats its behavior at regular intervals dependent on the 
statistical parameter. The interval time is the longest for 
fermions and is the shortest for bosons. The occupation 
distributions of natural orbitals do not change qualita- 
tively with the time evolution but the occupation of the 
lowest natural orbital evolves with time for anyon gas 
deviating the Fermi limit. When the harmonic trap is 
turned off it increases and gradually arrives at the biggest 
value. When the harmonic trap becomes weaker rather 
than being turned off, it shall fluctuate with the time 
evolution. 
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